Large nonlinear absorption and refraction coefficients of carbon nanotubes estimated 

from femtosecond z-scan measurements 



O 

o 

(N 
O 

Q 



o 

CZ3 



N. Kamaraju, Sunil Kumar, and A. K. Sood* 

Department of Physics and CULA, Indian Institute of Science, Bangalore 560012, India 

Shekhar Guha 

Materials Directorate, WPAFB, Dayton, OH 45433, USA 

Srinivasan Krishnamurthy 
SRI International, Menlo Park, CA 94025, USA 

C. N. R. Rao 

CPMU, Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore 560064, India 

(Dated: February 2, 2008) 

Nonlinear transmission of 80 and 140 femtosecond pulsed light with 0.79 /im wavelength through 
single walled carbon nanotubes suspended in water containing sodium dodecyl sulphate is studied. 
Pulse-width independent saturation absorption and negative cubic nonlinearity are observed, re- 
spectively, in open and closed aperture z-scan experiments. The theoretical expressions derived to 
analyze the z-dependent transmission in the saturable limit require two photon absorption coefficient 
Po ~ 1.4 cm/MW and a nonlinear index 7 ~ —5.5 x 10 -11 cm 2 /W to fit the data. 
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Single walled carbon nanotubes (SWNTs) have been 
studied for numerous applications including third or- 
der optical nonlinearity. These applications include na- 
noelectronics, gas and bio sensors, field emission dis- 
plays, saturable absorbers for passive optical regenera- 
tion, mode-locking and THz optical switching. The sus- 
pended carbon nanotubes show optical limiting due to 
nonlinear scattering, micro-plasma formation and subli- 
mation in the nanosecond regime In the femtosec- 
ond (fs) regime where heating does not play a role in 
the nonlinear transmission, an enormously large third- 
order susceptibility (Im\^ ~ 10 -6 esu) has been re- 
ported by resonantly exciting at the first inter sub-band 
energy levels (Su) of semiconducting SWNTs The 
origin of this large nonlinearity is assigned to the co- 
herence effect, rather than the incoherent or saturable 
absorption, because the measured nonlinearity decreased 
rapidly when the wavelength is changed away from the 
first band gap£ This interpretation for the large nonlin- 
earity is being debated*^ Other experiment a 10 ' 11 using fs 
pulses have reported a mixed variation with wavelength- 
much smaller valued for Imx'- 3 -' ~ 10~ 10 esu at 1.55 \im 
and a fairly large valued of ~ 1CP 7 esu at 1.33 \im. The 
origin of the nonlinearity had not been identified. In this 
letter we report the results of closed aperture (CA) and 
open aperture (OA) z-scan measurements carried out in 
suspension of almost isolated SWNTs at a wavelength of 
0.79 fim (energy hui of 1.57 eV). The OA z-scan shows 
saturated absorption and CA z-scan reveals negative cu- 
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bic nonlinear refraction. We have developed a theoreti- 
cal model incorporating saturated absorption along with 
nonlinear absorption and refraction to derive the trans- 
mission in both the OA and CA z-scans. The theoretical 
analysis of our results obtained with two different pulses 
with full-width at half-maximum (FWHM) of 80 fs and 
140 fs clearly identifies two photon absorption (TPA) as 
the source of nonlinearity and the TPA coefficient (3 is 
~ 1.4 cm/MW. This translates to a fairly large value of 
~1 x 10~ 9 esu for Imx^ 3 ' nonlinearity. 

A dispersion of SWNTs (0.4 mg) and 1% of sodium do- 
decyl sulphate (SDS) in water (1 ml) was sonicated for 
5 hours and the resultant solution was found to be well 
dispersed. This solution in 1 mm path length quartz cell 
was used in our experiments. The SWNT sample used 
in our experiments is pristine and contains two diame- 
ter distributions at 1.41 nm and 1.58 nm as confirmed 
by the radial breathing Raman modes at 160 cm -1 and 
177.7 cm -1 , shown in the inset of Fig. 1. 

Nonlinear transmission studies were carried out at 
1.57 eV using Ti: Sapphire Regenerative femtosecond 
amplifier (Spitfire, Spectra Physics). The chosen pho- 
ton energy is nearly resonant with the first interband 
transition energy in metallic tubes, M\\ and off- resonant 
with the second interband transition energy of the semi- 
conducting tubes Si2- From the absorption spectra of 
the dispersed nanotubes shown in Fig. 1, we infer that 
the absorption coefficient at this wavelength is about 
5.6 x 10 4 cm" 1 . The FWHM pulse width of the am- 
plifier output was 50 fs at a repetition rate of 1 kHz. 
Near the sample point, the pulse from the amplifier was 
found to be broadened to 80 fs. For the experiments 
done with 140 fs, we stretched pulse by adjusting the 
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FIG. 1: Optical absorption spectrum of SDS suspended Sin- 
gle walled carbon nanotubes in water. £jf 2 (-Eu) are the 
second (first) interband transition energy of semiconducting 
(metallic) nanotubes. The inset is the Raman spectrum show- 
ing radial breathing modes 



compressor of the amplifier. Two Si-PIN diodes (one for 
the signal (B) and the other for reference (A)) triggered 
at the electronic clock output (1 kHz) from the ampli- 
fier are used for the data acquisition and the difference 
between B and A was collected using a lock-in amplifier 
(SRS 830), averaged over 300 shots. This difference data 
was then converted into actual B/A signal in a personal 
computer. The SWNT dispersion in 1 mm thick cuvette 
was translated using a motorized translation stage (XPS 
Motion controller, Newport) over the focal region. The 
intensity of input beam was varied from 150 MW/cm 2 
to 6.2 GW/cm 2 . In the OA z-scan, all the light was col- 
lected by using a collection lens in front of the diode. The 
measured (and normalized) transmission data as shown 
in Fig. 2(a) clearly demonstrates the saturable absorp- 
tion where the transmission is enhanced at focus (z = 0). 
For CA z-scan, we kept an aperture of diameter 3.6 mm 
in front of the diode B and the measured transmission 
(and normalized) data as shown in Fig. 2(b) indicates 
photo-carrier induced reduction in the refractive index. 

The procedure to calculate the transmission in z- 
scan experiments has been well describe d 12 ' 13 ' 14 ! 15 in 
the literature. For the optical limiting case, fairly ac- 
curate solutions have been used to explain the z-scan 
resultsi 12 ' 13 i 14 ' 15 However, for the saturation absorption 
case, the solutions are obtained either qualitatively^or 
in the limit where the intensity is far smaller than the 
saturation intensity^ To quantify the basic mechanisms 
responsible for the nonlinear absorption and refraction, 
we accurately solve the rate equation, but modified^ for 
saturation absorption: 
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where evo is the one photon absorption (including intrinsic 
free carrier absorption) coefficient, /3q, is the fundamental 



TPA coefficient, and cto is free carrier absorption (FCA) 
cross section. I s is the parameter that characterizes the 
saturation absorption. The intensity (I) at the radial 
position r, time t, the position £ in the sample, and the 
location of the sample z is denoted as I(z, r, t). The 
generated photo-carrier density AN depends on both ao 
and (3q. In our experiment, the maximum value of / is 
~ 6 GW/cm 2 and ao is = 5.6 x 10 4 cm -1 . Even with 
a large value for /3 ~ 10~ 6 cm/W, the dominant source 
of carrier generation is one photon absorption. Since the 
carrier decay time is much longer than the pulse width 
tq, Eq. 1 takes the form: 
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The boundary condition required to solve Eq. 2 is the 
input intensity which is assumed to be a Gaussian and 



l(z,0,r,t) =I 



w 



w(z) 



exp 



2r 2 
w(z)" 2 



exp 



t 2 



To- 



il) 



wq is the beam waist at the focus, w(z) = 
woy/l + (z/zq) 2 is the beam radius at z, zq — ttwq/X 
is the diffraction length of the beam, To is the half 
width at e~ 4 of the maximum of the pulse, and A is 
the wavelength. The intensity at the exit side of sam- 
ple, I(z, L, r, t) is obtained from the analytical solution 19 
to Eq. 2 and integrated over all r and t to calculate the 
transmitted energy. The transmission in the OA z-scan 
experiment, Toa{z) is simply the ratio of transmitted 
energy to the incident energy. The solution to Eq. 2 de- 
pends on two parameters e g and I s , which can be varied 
to fit the measured OA z-scan data. Since several sets of 
/3 e ff and I s can fit the data, we will choose appropriate 
set that fits CA data as well. 

The transmission for CA case is more complicated as it 
requires the phase of electric field (E) in addition to the 
intensity at the exit surface. 12 ' 14 ! 15 The phase at £ — L 
is different from that at £ = because of the change in 
refractive index, An, caused by light absorption. To a 
first order in J, the An is simply 7/, where 7 is the cubic 
nonlinear index. 

The cubic nonlinear refraction coefficient, 7, in general 
has contribution from the photo-generated carriers, tem- 
perature change, and bound electrons. Following the well 
established procedure^, we calculate the electric field in- 
side the sample. For thin samples where photo-carrier 
generation is uniform, the change in phase from inci- 
dent to exit surface is Acj) = kj I(()d(, where k is 
the wavenumber and /(£) is the solution of Eq. 2. For 
a given r and t, the intensity at any £ inside the sam- 
ple is obtained analytically, fitted to polynomial series in 
I(z, 0, r, t), and substituted for A<f) to get, 
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The electric field at the exit surface is then, 



E(z, L, r, t) = v / I(z,L,0,0)e 



r T i ikr 2 
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where R(z) = z[l + z 2 jz 2 \ is the radius of curvature of 
the wave- front at z. As before^, e lA ^ is expanded in infi- 
nite series and Gaussian decomposition method^ 3 - is used 
to obtain the field pattern at the aperture which is at a 
distance d away from the exit surface. In our calculations 
we found that the convergence is achieved with first four 
terms in the expansion. We get normalized z-scan trans- 
mittance, Tca(z) by integrating \E(z, L + d, r, t)\ 2 over 
all t from — oo to +oo, and r from to aperture radius 

r a , then dividing it by ir?T w 2 I [l - exp( = -i^)]/2. 
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FIG. 2: Normalized transmittance data (open circles) in 
(a)OA z-scan and (b) CA z-scan Ql-ea;p(^£*)]=0.72). The- 
oretical fit (solid line) is obtained with (3 e g — 1.4 cm/MW , 
I s = 30 GW/cm 2 (for OA) and 37 GW/cm 2 (in CA) 



From the difference between the normalized peak 
and valley transmittance in CA z-scan experiment and 
aperture's linear transmittance, value of 7 can be 
calculated^ 2 -. From our measured CA z-scan experi- 
ments (shown in Fig. 2(b)), we obtain a value of 
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—5.5 x 10 -11 cm /W for 7. A simple calculation^ using 
the value of —3 x 10~ 21 cto~ 3 for (dn/dN) predicted 21 
in wide bandgap semiconductors also yields a value 
—3 x 10~ n cm? /W for 7. However our measured value 
for 7 is about two orders less compared to the value pre- 
dicted by Margulies et al for SWNTs^ 2 -. Using a value 
of —5.5 x 10 -11 cm? jW for 7, we have then varied /3 e ff 
in the calculation of Tca(z) to fit the data. We found 
that both OA and CA data can be fitted simultaneously 
with one set of parameters as shown by solid line in Fig. 
2. With larger values of /3 e ff, a good fit to CA data 
could not be achieved for any value of I s (the same is 
true for OA data too). Although good fit is possible 
with much smaller [3 e g, the I s required to fit CA differs 
considerably (by more than an order of magnitude) from 
that needed to fit OA data. A consistent set of parame- 
ters for 7, [3 e ff and I s respectively, —5.5 x 10~ n cm 2 /W, 
1.4 x 1CT 6 cm/W, and 30 GW/cm 2 , fit both OA and 
CA data well as shown by solid lines in Fig. 2. As seen 
from Eq. 3, the /3 e ff has two contributions -TPA and 
FCA. Noting that the FCA contribution depends on the 
pulse width, we repeated both CA and OA z-scan with 
FWHM width of 140 fs to evaluate the relative strength 
of these two contributing mechanisms. We found that 
z-scan data with 140 fs pulse width is identical to that 
obtained with 80/s pulse width. This clearly indicates 
that FCA cross section is extremely small and the funda- 
mental TPA (/3o) is the dominant mechanism for nonlin- 
ear absorption. The predicted value of lAcm/MW for 
Po in CNTs is two-to-three orders of magnitude larger 
than that in wide band gap semiconductors at 1 to 2 fim 
wavelength. This large value of nonlinearity (both (3q 
and 7) arises mainly as the consequence of the one di- 
mensional motion of delocalized-electron cloud along the 
nanotube axis . 23 ' 24 Although the large value of /?o at the 
wavelength, where «o is also large, makes it less useful for 
optical limiting applications, the underlying mechanism 
responsible for /3q could be operative even at the forbid- 
den wavelength. With our increasing ability to tune the 
band gap with nanotubes' radius, CNTs offer an interest- 
ing possibility for enhanced nonlinearity in near visible 
to short wave infrared wavelength region. 
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